arXiv:1509.02154vl [cs.RO] 7 Sep 2015 


Imperial College 
London 


Imperial College London 
Department oe Computing 


Self-calibration of a differential wheeled robot 
using only a gyroscope and a distance sensor 


Author: Supervisor: 

Carlos Garcia-Saura Prof. Andrew J. Davison 


Submitted in partial fulfillment of the requirements for the MSc degree in 
Computing Science (Specialism in A.I.) of Imperial College London 


September 2015 





Abstract 


Research in mobile robotics often demands platforms that have an adequate balance 
between cost and reliability. In the case of terrestrial robots, one of the available op¬ 
tions is the GNBot, an open-hardware project intended for the evaluation of swarm 
search strategies. The lack of basic odometry sensors such as wheel encoders had 
so far difficulted the implementation of an accurate high-level controller in this plat¬ 
form. Thus, the aim of this thesis is to improve motion control in the GNBot by 
incorporating a gyroscope whilst maintaining the requisite of no wheel encoders. 
Among the problems that have been tackled are: accurate in-place rotations, mini¬ 
mal drift during linear motions, and arc-performing functionality. Additionally, the 
resulting controller is calibrated autonomously by using both the gyroscope module 
and the infrared rangefinder on board each robot, greatly simplifying the calibration 
of large swarms. The report first explains the design decisions that were made in 
order to implement the self-calibration routine, and then evaluates the performance 
of the new motion controller by means of off-line video tracking. The motion accu¬ 
racy of the new controller is also compared with the previously existing solution in 
an odor search experiment. 

Keywords: Mobile robotics, dijferential wheeled robot, motion control, self-calibration, 
continuous-rotation servomotor, PID auto-tuning, multimodal sensing, MEMS gyro¬ 
scope, IR distance sensor, open-source 
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Chapter 1 


Introduction 


The starting point of this project is the GNBot swarm robot platform that was pre¬ 
sented in “Design principles for cooperative robots with uncertainty-aware and resource- 
wise adaptive behavior”^, and “Cooperative strategies for the detection and localiza¬ 
tion of odorants with robots and artificial noses’’^. Figure [T7T] describes the most 
relevant characteristics of these robots. 



Figure 1.1: One of the robots from the GNBot swarm used for the project. The electronics 
are based on the Arduino MEGA board. A custom shield contains all of the multimodal 
sensors, which provide the robot with gas sensing capabilities (TGS-2900), as well as dis¬ 
tance sensing (GP2Y0A21YK0F), temperature & humidity (DHTll), and light intensity 
sensing. This project additionally incorporated a gyroscope module (MPU6050). 
The main actuators are two continuous-rotation servomotors, and each robot also has 
a multicolor RGB LED and a piezoelectric speaker. Einally the wireless communication 
layer is based on ZigBee, and the green/yehow top markers allow for external video 
tracking of the swarm. From IfTHSH . 

The GNBot is a differential wheeled robot, which means that it has two active wheels 

sharing the same axis, and each of them is actuated independently. There is also a 
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passive, free turning ball caster that acts as a third stand. If the desired trajectory 
is a straight path, the velocities are set to equal magnitude; if a rotation is required 
instead, then different velocities can be applied to each motor. 

Using this setup it is possible to estimate the motion of the robot in the XY plane 
given its starting position (Pxya), the rotational velocity of each motor (uil, ujr), the 
diameter of the wheels Dr) and their separation (d). This estimation process 
is known as dead reckoning. 

The accuracy of dead reckoning relies on how dependable the measurements used 
in the calculations are. These include wheel separation and their diameter, but most 
importantly the actual rotational velocities of each motor. For instance, small devia¬ 
tions from the theoretical velocities can turn a desirably linear trajectory into an arc. 
An example of this issue is displayed in Figure 11.21 



-150 -100 -50 0 50 

Position in X [cm] 


Figure 1.2: Example of an odor search strategy based on Levy walks, as performed by 
the original robot controller. The search is defined at a high level as random in-place 
rotations combined with linear trajectories whose lengths vary according to a heavy¬ 
tailed probabilistic distribution. It can be appreciated that the segments that were 
supposed to be linear, are actually performed as arcs (there is drift in yaw). In 
this experiment, the target odor source (red marker) took 12 minutes to be found. The 
ground truth trajectory was recorded with a webcam, post perspective correction and 
color marker tracking with the OpenCV software library. 


^https://en.wikipedia.org/wiki/Dead_reckoning 
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Chapter 1. Introduction 


This thesis has tackled the improvement of those low-level control routines in order 
to minimize motion drift, as well as the implementation of a distance-based abstrac¬ 
tion layer that facilitates the use of the GNBot platform in practical applications. It 
also explores the automatic calibration of the robots using only on-board sensors (a 
gyroscope and a distance sensor). 

The following two sections are an overview of the project’s approach towards the 
improvement the motion controller; Chapter [2] provides an in-depth explanation of 
the self-calibration algorithm; Next, Chapter 0 evaluates the resulting controller; 
Finally, Chapter [4] analyses the outcome of the project and suggests future research 
paths. 


1.1 Background and motivation 


Undesired drift in the estimation of position, very characteristic when using dead 
reckoning in mobile robotics, can be reduced with the incorporation of basic odom- 
etfyd such as rotational encoder^. Wheel encoders are a very common choice given 
their simplicity and reduced cost, but they do have some drawbacks: 

• Incorporating wheel encoders into a previously existing design is often nontriv¬ 
ial, as the process involves hardware modifications near critical moving parts. 

• Electrical, optical or magnetic wheel encoders can be susceptible to dirt; plac¬ 
ing them near the wheels may make periodic maintenance necessary. 

• Most importantly, wheel encoders cannot easily detect whether there is any 
slipping between the wheel and the ground surface. 

For these reasons the incorporation of wheel encoders in the CNBot has been dis¬ 
missed. Instead, the first trials towards the reduction of yaw drift included the use of 
an electronic compass. The goal was to obtain a global, unbiased measure of head¬ 
ing (yaw) by measuring the earth’s magnetic field. Unfortunately the compass was 
deemed unusable for indoor environments during the first trials of the CNBot 11211. 
since the magnetometer measurements are substantially distorted by nearby metals 
(i.e. buried cables) commonly present indoors. 

Another option for odometry is the gyroscope, which is the approach studied in this 
project (see Fig. II.3D . Electronic gyroscopes are a form of inertial sensor that can 
measure relative rotations in any orientation, without requiring interactions with 
external moving parts. As opposed to wheel encoders, the physical incorporation of 
a gyroscope into the CNBot would be quite straightforward, and would also have 
the advantage of being capable of easily detecting wheel slipping. 


^Odometry: Use of sensory data to improve motion estimates 

^Rotational encoder: A sensor that can measure wheel rotations accurately and in real time 
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Figure 1.3; Detail of the gyroscope module installed in the GNBot. The selected board 
is based on the MPU6050 module, which is Arduino-compatible. The module fits on 
the Pc socket of the GNBot without any additional modification. An orange wire was 
soldered to the interrupt pin, but it was not needed in the end, as polling was used to 
read values from the module. 


The use of inertial sensors for position estimation has been a matter of research for 
many decades. The last few years in particular have seen the development of smart¬ 
phones, tablets, intelligent cameras and other devices in need of inertial sensors, 
which has greatly boosted the industry. The technology has been able to achieve high 
miniaturization, high accuracy and low power consumption, whilst substantially re¬ 
ducing costs, yielding the excellent inertial sensors that are now widely available in 
the market. 

In robotics. Inertial Measurement Units (IMUs) are quite often employed for navi¬ 
gation in different environments. For instance, the combination of IMUs with GPS 
allows for a remarkably accurate 3D localization in aerial robotics [|3l] . These tech¬ 
niques have also been demonstrated in ground robots without GPS||4ll5D, as well as 
in underwater robots without encoders Il6l [71 [8t]. More recently, 3D position track¬ 
ing with inertial sensors has also been employed to correct the effect of the rolling 
shutter of RGB and RGB-D iTTOfI camera sensors. 

Many authors use Kalman filter^ in their motion controllers. Kalman filtering is 
a probabilistic method for sensor fusion, which means it can be used to efficiently 
combine measurements from the inertial sensor with other forms of odometry such 
as wheel encoders, a compass, or even barometers; this way it is possible to achieve 
very robust motion controllers with minimal drift ||4l [2, [TTl [T2ll . In fact, most of 
the commercially-available IMUs, such as the one selected for this project, already 
contain internal Kalman filters that can integrate measurements from a gyroscope, 
an accelerometer or a compass. 

^https://en.wikipedia.org/wiki/Kalman_filter 
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Chapter 1. Introduction 


For the final evaluation of the positioning accuracy of the controller, the project took 
advice from work by J. Borenstein et al. IHTI [TBII , where the robots were programmed 
to describe square trajectories that were first logged and then compared against 
ground truth. 


1.2 Contribution overview 


At the start of this project, the motion control of the GNBot robots was speed-based 
and open-loop (with no wheel encoders). This resulted in very inaccurate motions 
that drifted rapidly from the desired path. For example, when a robot was com¬ 
manded to describe a straight trajectory, an arc path would be observed instead. 
In-place rotations were also quite inaccurate. 

Thus the goal of this project has been to improve the motion accuracy of these robots 
and to provide a high-level distance-based controller. In particular the following 
matters have been tackled: 


• Incorporation of an inertial sensor (MPU6050) into the GNBot design, as well 
as the implementation of low-level functionality that reads yaw, pitch and roll 
from the gyroscope module with minimal drift. 

• Analysis of the servomotor response curves. Since there are no wheel encoders, 
the rotations were tracked using the gyroscope, by independently rotating each 
wheel over the entire range of velocities. 

• Implementation of a PID controller for accurate yaw regulation; Additionally 
the constants are auto-calibrated with a process based on the Ziegler-Nichols 
method. 

• Analysis of the infrared sensor response curve, and linearisation using expo¬ 
nential curve fitting. 

• Analysis of the nonlinear effect of the PID yaw controller over the robot’s linear 
velocity for the whole range. 

• The results of these analysis have been integrated into a high-level distance- 
based motion controller that supports linear motions and arcs. Additionally, 
the calibration of each robot is done autonomously by using only on-board 
sensors (the gyroscope and the IR rangefinder) and requiring minimal user 
interaction. 

• Implementation of an off-line video tracking method capable of reliably record¬ 
ing ground truth robot paths. The process involved perspective & distance 
correction using the OpenCV library. 
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• Evaluation of the new controller by making the robots perform lines, circles 
and squares of known dimensions; Also the system was compared against the 
previously existing solution in an odor search task based on a wall-bounce 
strategy. 


The new motion controller has a linear drift in the order of 6cm/m in the low velocity 
settings (2 to 6cm/s) and around 9cm/m at higher velocities (6 to lOcm/s). The 
rotational drift is in the order of 2deg/min when the robot is moving. 

In summary, this project has provided the high-level functionality needed in order 
to achieve accurate control over the GNBot robots. The self-calibrating nature of the 
approach also facilitates its use in large robot swarms. 

The outcome of the project is open-source (Attribution-ShareAlike 4.0 International!) 
and can be accessed in the following GitHub respository: 


https://github.com/carlosgs/GNBot 


^http://Greativecommons.org/licenses/by-sa/4.0/ 


6 






Chapter 2 


The self-calibration algorithm 


The original controller on-board the GNBot was velocity-based and its input units 
were arbitrary - it did not have an internal mapping between the motor’s inputs and 
real-world units (i.e. cm/s). It was also lacking methods for basic calibration, in 
part due to the absence of any form of odometry. Altogether, these facts lead to large 
amounts of drift. Far from being a constant that could have been compensated by 
basic calibration, drift varied with the velocity of the robots -not only in magnitude 
but also in direction-, so straightforward calibration was not an option (see Figure 
I2.1D . This project has tackled the problem by incorporating a gyroscope for odome¬ 
try. The advantages of this approach in contrast with the use of wheel encoders were 
already discussed in Section [TTTl 

The selected Inertial Measurement Unit (IMU) is the MPU605cE] chip, which is based 
on MEMS technology (Micro-Electro-Mechanical Systems). It is commonly available 
as an Arduino-compatible break-out board that conveniently fits on the PC socket 
of the GNBot without any additional modifications (c.f. Figure 11.31) . 

The gyroscope module has the same PC pinout as the magnetometer that was orig¬ 
inally present in the robot, but the routines that handle both modules are very dif¬ 
ferent. In fact, the example code for the inertial sensor (a part of the PC Device 
Librar^ depends on the use of hardware interruption^ rather than polling. The 
use of interruptions was highly undesirable since the main processor needs to han¬ 
dle time sensitive tasks such as sensing and radio communications as well as motion 
control, so the gyroscope library was modified to use a polling scheme. 

The MPU6050 incorporates a Digital Motion Processor (DMP) that can be used to 
offload computations from the robot’s main processor. These computations include 

'http://www.invensense.coni/products/motion-tracking/6-axis/mpu-6050/ 

^http://www.i2cdevlib.com/devices/mpu6050 

^Interrupt: Signal to the processor indicating that some event needs immediate attention 

'^Polling: As opposed to interruptions, polling actively samples the status of an external device 
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the task of filtering and integrating the raw values from the on-chip acceleration 
sensors, in order to output ready to use yaw, pitch and roll. In the GNBot controller 
the raw values are sampled and processed at 200Hz by the DMP co-processor. 

The gyroscope module also has an internal FIFcH buffer that allows for reliable high- 
frequency data readings (up to 200Hz), though in this project it was deemed unnec¬ 
essary since the sample rate is low (< 50Hz). Once the IMU had been incorporated 
into the GNBot, it was possible to read the yaw, pitch and roll of the robot back into 
the computer (see Figure IZTl) . 


Observed trajectory for distinct linear velocities 



Deviation from a straight path (yaw), measured with the on-board gyroscope 



Figure 2.1: Initial experiments to analyse drift in the original controller. The robot was 
commanded to describe a straight path by setting both motors to the same velocity 
during 3 seconds. The trajectory (upper panel) was recorded with a ceiling camera 
for three different velocities, while yaw (lower panel) was being logged using the on¬ 
board gyroscope. It can be observed that drift is accurately tracked by the empirical 
yaw measurements. The input speed setting is arbitrary since both continuous-rotation 
servomotors are commanded with the built-in Servo.write{[degrees]) functions from the 
Arduino IDE, which do not have a direct mapping with real world velocities. 

At this point, a basic proportional yaw controller was implemented in order to 
^FIFO: “First In First Out” queuing policy 
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Chapter 2. The self-calibration algorithm 


demonstrate the viability of the system. The controller could now successfully cor¬ 
rect the direction of the robot so it followed a straight path, even after external 
perturbations had been applied (i.e. rotating the robot or placing obstacles). 


With those results as a motivation, the design of a high-level self-calibrating motion 
controller was tackled. An overview of the resulting calibration routine is shown in 
Figure [2^- Further sections will describe the technical details in greater depth. 



Figure 2.2: Full self-calibration procedure, (a) The robot is first placed in front of a wall, 
and remains static while the gyro drift is compensated, (b) The minimum motion thresh¬ 
old is found by gradually increasing the input of each motor until a rotation is perceived 
(c) The motors are then independently run at their maximum speed setting in order 
to record the corresponding real-world rotational velocity, (d) Next, the yaw PID con¬ 
troller is heuristically calibrated by analysing the oscillations for different parameters, 
(e) Finally, the robot uses the distance sensor to measure its approach velocity towards 
the wall, creating a map between the input values {rad/s) and the resulting real-world 
linear velocity (cm/s). 

Technical details can be found in sections: (a) 12.11 (6 & c) 12.21 (d) 12.31 (e) 12.41 
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Chapter 2. The self-calibration algorithm 


2.1 Minimization of gyroscope drift 


Most electronic gyroscopes calculate rotational magnitudes by integrating measure¬ 
ments from rotational accelerometers over time. Ideally this integration process 
would consistently yield the same outputs when a sequence of rotations is applied, 
but instead, rotational accelerometers often have inaccuracies (i.e. bias, sensitivity 
limits, variability with temperature, etc.) that cause drift in the integration process. 
This would imply for instance that the yaw measurement of a static robot would 
erroneously vary over time. 

Gyroscope drift is often minimized with either a proper calibration of the accelerom¬ 
eter bias, with signal filtering or by pausing the integration process while no motion 
is detected. Fortunately, the DMP co-processor of the MPU6050 already implements 
these features very efficiently (see Fig. 12.3D . 


In order to measure yaw drift rate (whose units are degrees per minute) it was nec¬ 
essary to differentiate the yaw measurements provided by the gyroscope (radian 
units). This was achieved by accurately timing distinct yaw measurements with the 
function millisQ from ArduincH. 


driftjrate 


yaw{t2) - yawifi) 
^2 ~ 


ISOdegf 

Tirad 


Imin 

- deq min 

60s 


( 2 . 1 ) 



Figure 2.3: Analysis of gyroscope drift (with static robot). The gyroscope has an initial 
transient that must be respected, during which the internal DMP (Digital Motion Proces¬ 
sor) performs self-calibration to account for drift. The module is fully calibrated after 23 
seconds of the robot being static, when drift becomes inferior to 0.5 deg/min (step A in 
Fig. 12.2D . It must be noted though that the magnitude of this drift will not be constant, 
but is instead greatly affected by the robot’s motions. 


^https://www.arduino.cc/en/Reference/Millis 
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2.2 Calibration of the motors 


This section tackles the calibration of the continuous-rotation servomotors, which 
are actuators that generate a rotational velocity that is proportional to a PWM input 
signal (Pulse Width Modulation). The relationship between the input and output 
magnitudes is arbitrary, so a real-world mapping needs to be learned first. 

To better understand the calibration process it helps to be familiar with the following 
notation: 



Figure 2.4: Diagram of the robot’s top view, displaying the notation used. The grey arrow 
points towards the front of the robot. Dl and Du are the wheel diameters, and d is their 
separation. Vl and Vr are the contributions of each motor to the robot’s linear velocity, 
and a is the yaw measurement (positive for clockwise rotations). 

Most mobile robotic platforms have encoders that simplify the calibration process 
by measuring the real-world rotational velocities of the wheels. In the case of the 
GNBot there are no wheel encoders, so a different approach was needed. 

The implemented method can measure the response curve of each motor indepen¬ 
dently by using only the gyroscope. This process is described in Figure [231 



Figure 2.5: Method employed to measure lor, the contribution of the right motor to the 
robot’s rotational velocity uj. A constant velocity is first applied to the right wheel {Vr = 
C) while the left wheel is fixed {Vl = 0). This results in a rotation of the robot around 
the ground-contact point of the left wheel. The gyroscope is then used to measure the 
rotational velocity ujr = to. Afterwards the same process is repeated for the left motor. 
This method allows to separately evaluate the response curves of each motor without 
requiring wheel encoders. 
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Chapter 2. The self-calibration algorithm 


The angular velocity u, directly measured by the gyroscope, is coupled with the 
rotational speed of both motors with constants Kl & Kr, that account for Dl, Dr 
(wheel diameters), and d (distance between wheels). Using this technique it was 
possible to record the response curve of each motor by independently performing a 
sweep over the full range of velocity inputs. The result is shown in Figure [Z6l 


Velocity response curves for distinct wheel diameters 

Dl = 75mm, Du = 65mm Di = Dr = 65mm Di = 65mm, Dr = 75mm 



Figure 2.6: Measured motor response curves. The panels show the response curves of the 
left (red) and right (blue) motors for different wheel sizes. The input of the continuous- 
rotation servomotors is PWM (Pulse Width Modulation), and the output (w) is the ro¬ 
tational rate of the robot around its vertical axis (the setup is described in Fig. 12.5ft . 
Each trial was run for a single motor at a time, and the sweep took around 30 seconds 
per motor. Three regions can be appreciated: a flat dead-zone (uj = Orad/s), a linear 
response region (0 < w < 1.5rad/s), and saturation (w > 1.5rad/s, which has been 
highlighted in red). For the same wheel size -middle panel- the L/R response curves 
are almost identical. Robots with larger left or right wheels -left & right panels- instead 
have visibly different L/R response curves. The calibration routine implemented in this 
project accounts for these differences. 

Self-calibration of these curves has been implemented by means of linear fitting 
between the minimum and maximum velocities within the linear region. In first 
place, the robot measures the dead-zone of a motor (the minimum motion threshold) 
by gradually increasing its input velocity until a rotation is perceived. Then, the 
maximum velocity is measured by running the motor at the maximum speed setting 
within the linear range (PWMsilSOO or PWMsilbSO). These measurements are 
performed for the forward and backward directions of each motor (steps B & C in 
Fig.Q. 

The left and right motors contribute to the robot’s global angular velocity oj in an 
additive manner (oo = ool + ujr, i.e. setting ool = 0.5rad/s and ujr = 0.5rad/s 
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Chapter 2. The self-calibration algorithm 


yields u = 1. Or ad/s). Knowing this fact, and using the data from Figure it 
was possible to determine the input region that corresponds to a linear motion. 
This is when both rotational contributions cancel out so the robot does not rotate: 
u = ul + Ur = Orad/s ur = —ur. This region is represented in Figure [2?7l 
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Figure 2.7: Region of motor input values that correspond to a linear trajectory. The red 
and green data points are from robots with larger left or right wheels respectively, while 
the blue data points are from a robot with identical wheels. It can be observed that each 
of these curves has a different steepness: a greater slope indicates that the right motor 
needs to rotate faster in order to account for a larger left wheel, and vice-versa. 
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2.3 PID auto-tuning using the Ziegler-Nichols method 


At this point a basic motor controller had been implemented, so the rotational rate 
{rad/s) of the robot could be accurately commanded. Closed-loop yaw control was 
now possible with the implementation of a PID regulator (Proportional-Integral- 
Derivative) that combines the motor controller with real-time gyroscope measure¬ 
ments (Eqn. 2.2). The error measure is calculated as e{f) = yawtarget — yaw{t). 


output_rotationjrate{t) = Kpe{t) + Ki / e{T)dT + Kd-y-e [rad/s] 

In at 


( 2 . 2 ) 
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Chapter 2. The self-calibration algorithm 


A PID controller is based on user-tunable gains (Kp, Ki & Ka) that, when properly 
adjusted, can minimize the overshoot of the transient response; in the case of the 
GNBot, the PID is the responsible of achieving fast and accurate changes in yaw. This 
is represented in Figure 12.81 




Time [seconds] 

Figure 2.8: Response of the yaw controller to a step function (upper panel) and to a 
sinusoid (lower panel). The controller was commanded with time-varying yaw targets 
while the transitions were recorded using the gyroscope. It can be appreciated that a 
proportional controller either is too slow (cyan curve), has over-shoot (blue curve) or 
even ripple (red curve). The auto-tuned PID (green curve) does have some overshoot as 
well, but it is a much better approximation to the ideal response. 

Rather than manually specifying a fixed parameter set for the PID, this project tack¬ 
led the self-calibration of the controller with the popular technique proposed by 
Ziegler and Nichols [fT4ll . This heuristic tuning method is performed by setting the 
three gains (Kp, Ki & Kf) to zero, and then increasing the proportional term Kp 
until it reaches the ultimate gain at which the controller presents self-sustained 
oscillations of period T^. Ziegler and Nichols (ZN) then provide with the mapping 
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Chapter 2. The self-calibration algorithm 


between Ku & Tu and the three gains of the PID controller. In this project the mod¬ 
ified tuning values proposed by Tyreus and Luvben lfTSlI were used instead as they 
increased the robustness of the controller. Table 2.1 compares the gain specifications 
of both tuning methods. 



Kp 


to 

Ziegler-Nichols (ZN) 

Kjl.f 

Tu/2 

Tu/S 

Tyreus-Luyben (TLC) 

KuI2.2 

2.2Tu 

Tu/6.3 


Table 2.1: Prescribed PID gain values for the Ziegler-Nichols (ZN) & Tyreus-Luyben (TLC) 
heuristic tuning methods. Compared to the original values proposed by Ziegler & Nichols, 
the TLC rules tend to reduce the oscillatory effects and improve the robustness of the 
yaw controller. The tj and to values were converted into the standard PID form using 
Ki = Kp/Ti and Kd = KpTd- 


In order to provide a consistent PID calibration, the parameters Ku -ultimate gain- 
and Tu -oscillation period- needed to be measured accurately and reliably. An auto¬ 
mated method was implemented for this purpose (see Figure [2?9l) . 


PID auto-tuning i:)roccss: Empirical measure of the values and 



Figure 2.9: PID auto-tuning process: Empirical measure of the values Ku & Tu- The 
tuning is done as follows: In first place, Ki & Kd are set to zero, and Kp is set to a very 
large value that ensures oscillation (the value is such that the maximum motor speed is 
applied at 1 degree of error). Next, the proportional controller is perturbed by setting a 
45 degree angle as the target yaw. Then the algorithm checks whether the oscillation is 
self-sustained, in which case it updates Kp = Kp/l.b. This process is repeated until the 
oscillations become attenuated; at that point Ku & are finally registered (step D in 
Fig. [121). 
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Chapter 2. The self-calibration algorithm 


2.4 Mapping rotational velocities with linear motions 


Using the basic motion controller it was now possible to command the robot to 
perform linear trajectories quite accurately; but the inputs units were still rotational 
magnitudes related to yaw (ujl & ujr [rad/ s]). The next step was to create a mapping 
between these inputs and the real-world linear velocities of the robot (i.e. cm/s). 

In first place it was necessary to have a method that could accurately measure the 
robot’s motions. External tracking systems (i.e. VICOI\0 or a ceiling camera) are 
among the best solutions available for this purpose, but they were deemed too so¬ 
phisticated to be incorporated into the basic calibration routine. These methods 
were employed for the final evaluation of the motion controller instead (Chapter [3]). 

Rather than requiring an external tracking system, the problem of measuring the 
robot’s real-world velocity has been tackled with the use of the distance sensor on¬ 
board the GNBot. The sensor is the Sharp GP2Y0A21YK0^mhared rangefinder, and 
it has the voltage response curve analysed in Figure 12.101 



Figure 2.10: Analysis of the response curve of the infrared distance sensor. The figure 
displays the output voltage of the Sharp GP2Y0A21YK0F rangefinder and its variability 
with the distance between the sensor and a wall. Manufacturer specifications rate this 
module for distance measurements in the range 10cm < d < 80cm, and indeed it can be 
appreciated that the response outside that region is either flat (d > 80cm) or nonlinear 
(d < 10cm, highlighted in red). 

The approach undertaken towards the calibration of linear motions was the use of 
the IR distance sensor on-board the GNBot to achieve accurate real-world velocity 
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^http://www.vicon.com/products/camera-systems 

®http://www.sharpsma.com/optoelectronics/sensors/distance-measuring-sensors/GP2Y0A21YK0F 
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Chapter 2. The self-calibration algorithm 


measurements. Theoretically, the robot would be driven towards a wall at a constant 
speed, and then the robot’s velocity could be calculated by differentiating distance 
measurements over time. 

In order to achieve a fair degree of resolution in these measurements, the output of 
the IR sensor needed to be calibrated first by means of linearisation. The method is 
introduced in Figure 12.111 


IR rangefinder response curve (II) 




Figure 2.11: Linearisation of the IR rangefinder sensor response. The left panel shows the 
response curves that were recorded for three different GP2Y0A21YK0F distance sensors 
(red, green and blue curves). These curves were fitted with an exponential decay (cyan 
curve) in the range of highest sensitivity (10cm < d < 40cm). The right panel evaluates 
the accuracy of the fitted model, showing a maximum theoretical deviation of Ri 3cm. 

The fitting function in Figure [2)TT] is an inverse function with two fitting parameters 
K=21.6V/cm & C=t).lQV that were calculated using two data points located in the 
limits of the highest sensitivity range {di = lt)cm & d 2 =^Dcm), using the following 
equations: 


vUd) = Kf + C K = 


d2 


K 

c = v{d2) - -J- 
(12 


(2.3) 


The calibration of the IR distance sensor is not a part of the self-calibration routine 
since it requires user interaction; instead, the fitting curve can be generalised to 
every robot that uses the GP2Y0A21YK0F sensor model. 


This newly-calibrated distance sensing functionality was verified to work within the 
theoretically-calculated tolerances for different wall surfaces. Having this infrastruc¬ 
ture in place, the problem of linear velocity calibration could finally be tackled. For 
this purpose, the GNBot was commanded to move towards a wall at different veloci¬ 
ties while yaw was corrected by the closed-loop PID controller. The real-world linear 
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Chapter 2. The self-calibration algorithm 


approach velocity (v cm/s) was then recorded for multiple input velocity settings 
(cuc rad/s), yielding the results shown in Figure [2T2l 


Linear velocity calibration for two different surfaces 
Surface type: Wooden table Surface type: Polyurethane foam 



Figure 2.12: Linear velocity calibration for two different ground surfaces. The panels show 
how the real-world velocity of the robot varies over its input range, for both backward 
and forward motions (red and blue curves) and for two different ground surfaces (left 
and right panels). The green and magenta dashed lines represent a linear fitting of the 
velocity response of each material. It can be appreciated that the speed of the robot 
over foam (right panel) is lower for the same input values; this is due to friction effects. 
Friction also causes an asymmetry in the velocity responses for backward and forward 
motions; this is because the robot is not perfectly symmetric. The nonlinear region is 
a cause of the PID operating the motors outside their linear range (c.f. Fig. 12.61) : such 
effect could be minimized with velocity-dependent PID gains. 

Self-calibration was then implemented by means of a linear fitting that measures 
the real-world velocity at the highest velocity setting within the linear region. This 
measurement is performed four times and then averaged (step E in Fig. 12.2D . being 
the last step of the self-calibration process. Each of the calibration parameters are 
then stored into the non-volatile EEPROM memory built into the Arduino, so the 
tuning process is not required every time the robots are powered on. 

Finally, high-level distance-based control was implemented by integrating over time 
the current velocity setting; this way the robot can stop whenever a specified dis¬ 
tance has been reached. This is a very simple method that does not account for 
transients in velocity. Arc functionality was implemented by interpolating the tar¬ 
get yaw throughout a distance-based motion; specifying different initial/final yaw 
angles effectively converts a linear segment into an arc. 
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Chapter 3 


Evaluation of the motion controller 


The purpose of the controller described in Chapter [2] is to provide accurate high- 
level motion control for the GNBot robots. This chapter evaluates the real-world 
performance of this controller. 

In first place, a series of experiments were designed in order to be able to measure 
the different motion drift characteristics. For instance, the robots would be com¬ 
manded to perform squares and circles of known dimensions while external video 
tracking provided unbiased measurements of the actual trajectories. The tracking 
process is described in Figures \3A] and [3721 



Figure 3.1: Perspective correction for accurate video tracking. The video tracking system 
was calibrated using perspective correction with a square of known dimensions. This 
allows for a 1:1 mapping between image units and real-world distances in the XY plane. 
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Chapter 3. Evaluation of the motion controller 



Figure 3.2: How color markers are tracked with the OpenCV library. In first place the 
reference points for perspective correction need to be specified by using the cursor (1). 
Next, the marker color is selected (2), and then a threshold is applied (3). The algorithm 
finally looks for the largest blob, which yields the coordinates of the robot (4). From ||2]] . 


Using this tracking system it was possible to record the trajectories of each robot with 
a very high resolution. Next sections describes the experiments that were performed. 


3.1 The effect of gyroscope drift 

Though gyroscope drift is always compensated upon the initialization of the GN- 
Bots, such compensation is far from ideal, and as a consequence the robot’s motions 
present yaw drift over time. This effect is analysed in Figure [3^31 



Figure 3.3: Cumulative yaw error due to gyroscope drift. In this experiment, the robot 
was commanded to repeat a linear trajectory multiple times (alternating the target yaw 
between 0° and 180°) for both forward motion (red path) and backward motion (green 
path) at a velocity setting of 5cm/s. Using this data, the rotational drift of the controller 
has been calculated to be in the order of 2 deg/min in a moving robot. 
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Chapter 3. Evaluation of the motion controller 


3.2 Analysis of linear positioning accuracy 

Next, the general performance of the distance-based controller was evaluated by 
first commanding the robots to perform squares and circles and then comparing the 
results against ground-truth dimensions. Square trajectories were implemented with 
straight line segments interleaved with in-place rotations of 90°. This functionality 
yielded the results shown in Figure 13.41 



Position in X [cm] 

Figure 3.4: Evaluation of the linear accuracy of the distance-based controller. A GNBot 
robot (lower-right corner) was commanded to describe square trajectories of different 
sizes while its path was recorded using the vision-based method described in Figures [STTI 
& 13.21 The red path corresponds to a purposely-miscalibrated robot (it was calibrated 
with a different set of wheels). The cyan and magenta paths correspond to robots with 
different wheel sizes that were properly self-calibrated. Finally, the green path corre¬ 
sponds to a self-calibrated robot with same wheel diameters. Linear drift has been cal¬ 
culated to be around 3cm/m. This was done by scaling up the maximum deviation from 
the lower-left corner of the largest square trajectory - about 1.5 cm for the calibrated 
robots. It can be appreciated that the shape of the squares is correctly maintained thanks 
to the gyroscope -which is always calibrated upon start up-, and the square dimensions 
only fall out of the specification for the miscalibrated robot (red path). 

In order to simplify the video tracking process, measurements were conducted in 
one single continuous experiment for each robot. In the case of square trajectories, 
the GNBot would perform the different square sizes one after another - without an 
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Chapter 3. Evaluation of the motion controller 


intermediate re-positioning of the robot. This generated an undesired drift in the 
starting position of each trial, which made manual alignment necessary. The post 
processing consisted in a affine translation from the starting points to the origin 
(x, ?/=0,0) and posterior re-setting of the initial yaw angle by means of a global 
rotation. This process effectively removed the effects of yaw drift between each trial 
(the effect had already been analysed in Figure [3^3]) . 

Circular motions were implemented with the arc functionality, by setting the initial 
and final arc angles to 0° and 360°. Path length was calculated as l=2Tir. This 
approach yielded the results shown in Figure 13.51 



Figure 3.5: Evaluation of the accuracy of arcs performed by the controller. A GNBot robot 
(lower-right corner) was commanded to describe circular trajectories of different sizes 
while its path was recorded. The cyan and magenta paths correspond to robots with 
different wheel sizes that were properly self-calibrated. The green path corresponds to 
a self-calibrated robot with same wheel diameters. Finally, the red path corresponds 
to a different robot that was also self-calibrated. The maximum drift from ground- 
truth diameters is around 8cm/m, though it can be appreciated that the circular shape is 
correctly maintained thanks to the gyroscope (which is always calibrated upon start-up). 
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Chapter 3. Evaluation of the motion controller 


3.3 How velocity affects the controller 


The implemented motion controller relies on the integration of velocities over time 
for its distance estimations. In Section 2.4 it was observed that there is a nonlinear 
region in the velocity response curves (c.f. Figure [27T21) . Since these responses are 
linearly fitted, the use of velocities outside the linear range will result in inaccurate 
distance estimations. The effect of velocity is analysed in Figures [3^ & 13 . 71 


Effect of velocity in the repeatability of square trajectories 



Position in X [cm] 



Position in X [cm] 

Larger right wheel 



-35 -30 -25 -20 -15 -10 -5 

Position in X [cm] 


Figure 3.6: Effect of velocity in the repeatability of square trajectories. The setup of these 
trials is similar to the one in Figure [3^ but in this case the size of the squares were set 
to a constant / = 30cm and performed at different speeds. The maximum drift is in the 
order of 6cm/m in the low velocity settings (2 to Qcm/s) and around 9cm/m at higher 
velocities (6 to lOcm/s). 
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Chapter 3. Evaluation of the motion controller 


Effect of velocity in the repeatability of circular trajectories 



Figure 3.7: Effect of velocity in the repeatability of circular trajectories. The setup of these 
trials is similar to the one in Figure [331 but in this case the circle diameters were set to a 
constant D = 30cm and performed at different speeds. The maximum drift in diameter is 
around 7cm/m in the low velocity settings (2 to 6cm/s) and around 8.5cm/m at higher 
velocities (6 to lOcm/s). 
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Chapter 3. Evaluation of the motion controller 


3.4 Performance in odor search experiments 


The GNBot robots were originally designed as a platform that could evaluate swarm 
search strategies in real-world environments [HI [2l] . The target of those searches are 
odor sources based on volatile compounds (such as ethanol) that can be detected by 
the gas sensor on-board each robot. 

Among the search strategies that have been so far evaluated with this platform are 
Levy-based search (c.f. Fig. II.2D and wall-bounce search. Unfortunately, at the time 
of those trials a calibrated motion controller was not yet present. This fact caused 
great uncertainty in the evaluation process: the calculated motion commands for a 
particular search were not being faithfully reproduced by the swarm. 

Some of those experiments have been re-run for this project - this time using the 
self-calibrated motion controller in order to evaluate the improvement. Trials were 
performed in the seminar room shown in Figure 13.81 



Figure 3.8: Room used for the odor search experiments. The odor search experiments 
required an effective area of approximately 16m^, an indoors environment free of dis¬ 
turbances and with uniform light conditions that facilitate video tracking. Chairs were 
moved to the back of the room, and the search arena was installed in the center (blue 
walls). Air conditioning was turned off during the experiments on order to minimize air 
currents. The robot swarm can be seen on top of the table, as well as the computer used 
to record tracking data. 

Wall-bounce search is a very primitive form of search where a robot is commanded 
to maintain a constant velocity until either an obstacle or the target are encountered; 
for obstacles, the robot must rotate away and continue with a different heading. The 
simplicity of the wall-bounce strategy makes it very sensible to the linearity of the 
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robot’s motions. For instance, a robot with enough yaw drift may rotate in circles 
over the same area without ever finding any obstacle that modifies its path. Given its 
sensitivity to drift, wall-bounce search is a very good candidate for the comparison 
of both motion controllers. The experimental setup is described in Figure [3T9l 



Square of known dimensions 
(for perspective correction) 


Odor sources (goal) 


Starting position 
' of the robots 


Figure 3.9: Setup of a wall-bounce search experiment (conducted with the new motion 
controller). The upper panel shows the starting robot positions, as well as the location of 
the odor targets. These are based on cotton pads impregnated with ethanol. The lower 
panel shows an overlay of the search path described by each robot. It can be appreciated 
that the trajectories are accurately performed as lines, which faithfully matches the high- 
level definition of the wall-bounce algorithm. 

The results of these search experiments are further discussed in the next section. 
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Chapter 4 


Conclusions 


This thesis has tackled the improvement of the motion controller on board the GNBot 
by incorporating gyroscope-based odometry that does not require wheel encoders. 
The outcome is a high-level controller that can faithfully perform in-place rotations, 
arcs, and linear trajectories. This controller is calibrated, so its inputs are real-world 
units (centimeters and degrees). Additionally this calibration process is automatic 
and uses only two on-board sensing capabilities: the gyroscope and the distance 
sensor. 

The general approach of the project has been to separately analyse the sources of 
motion uncertainty, and then try to find ways to minimise them (i.e. by observing 
a response curve and working only within its linear region). Many of the error 
sources could not be separated, so the calibration process had to progressively build 
a layered confidence base, which in the end allowed for complete self-calibration: 

• In first place, the output of the gyroscope (the MPU6050 inertial sensor) was 
analysed. It was found that the yaw measurement varied over time even when 
the robot was static, so this drift had to be compensated (Section [2.1D . 

• Then the response curves of the motors were analysed using the calibrated 
gyroscope. This was done by activating a single motor at a time and measuring 
the rotational velocities of the robot (c.f. Fig. 12.5D . The response curves were 
plotted for robots with different wheel diameters (c.f. Fig. 12.6D . These were 
then linearised, which allowed control over real-world units of rotational speed 
(Section [2^. 

• Using the calibrated motor functions, a PID yaw controller was implemented 
and tested in order to guarantee accurate yaw control (c.f. Fig. 12.8D . This 
allowed for fast and accurate yaw transitions. It also allowed for very accurate 
closed-loop linear trajectories (Section [2.3D . 

• The IR distance sensor on-board the GNBot (Sharp GP2Y0A21YK0F) was then 
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evaluated and linearised (c.f. Figs. 12.101 & [2.11D . Finally, the PID yaw con¬ 
troller was used in combination with the distance sensing capability in order 
to perform accurate velocity measurements. This way it was possible to mea¬ 
sure the linear velocity response curves of the robot (c.f. Fig. I2.12D . Basic 
fitting allowed for accurate velocity-based control with real-world input units 
(Section [Z4l) . 

The self-calibration process that has been implemented (c.f. Fig. 12.2D uses only the 
gyroscope and distance sensors on board the GNBot to automatically measure all the 
necessary parameters, without requiring any user interaction. The process simply 
requires manual placement of the robot in front of a wall, and takes approximately 
92 seconds. 

Video tracking was used for the evaluation of the controller’s accuracy (Chapter 
[3l) . Tracking was implemented using the OpenCV library, first performing perspec¬ 
tive correction, and then identifying the position of the color marker by applying 
a threshold (c.f. Figs 13.11 & 13.2D . This tracking method was used to compare the 
new controller with the previously existing solution, in a wall-bounce odor search 
experiment (see Figure 14.ID . 


“Wall bounce” search using the original controller (left) and the new self-calibration (right) 



0 100 200 300 400 500 600 700 800 

Position in X [cm] 


Figure 4.1: Comparison of a “Wall bounce” search using the original controller (left) and 
the new self-calibration (right). The figure shows the trajectories of four GNBot robots 
in two different odor search experiments. The starting positions are represented with 
triangle markers and the targets are represented as squares. It can be appreciated that 
the original controller (left panel) has large amounts of yaw drift. The new controller 
(right panel) faithfully matches the theoretical definition of wall-bounce search. 
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Chapter 4. Conclusions 


The tracking method also allowed to measure the effect of gyroscope drift (by hav¬ 
ing the robot repeatedly perform straight lines, as seen in Section IXTl) . the linear 
positioning accuracy (by performing squares and circles of different sizes, as seen in 
Section [321), and the effect of velocity in the controller (by performing squares and 
circles at different velocities, as seen in Section I3.3D . 

The new motion controller has a linear drift in the order of 6cm/m in the low velocity 
settings (2 to 6cm/s) and around 9cm/m at higher velocities (6 to lOcm/s). The 
rotational drift is in the order of 2deg/min when the robot is moving. 

In summary, this project has implemented the high-level functionality needed in 
order to achieve accurate motion control of the GNBot robots. The self-calibrating 
nature of the approach also facilitates its use in large robot swarms. The outcome of 
the project has been published as open-source in a GitHub repositor)0. 


4.1 Future work 


Using the new controller it will be possible to evaluate search strategies that de¬ 
pend on accurate position-based control (i.e. environment mapping, efficient area 
covering, flocking behaviours, etc). 

In regards to self-calibration, the studied method provides static tuning parame¬ 
ters that remain unchanged until a re-calibration is manually issued. Some authors 
have studied SLAM approaches that instead update the calibration estimates auto¬ 
matically in real time IfT^ [TTl . These techniques could be explored with the 
incorporation of other forms of odometry such as wheel encoders. 

The inertial sensor used in this project includes not only a gyroscope but also an 
accelerometer. It would be very interesting to incorporate these measurements into 
the motion controller. For instance, the accelerometer could be used in order to 
detect changes in ground surface and account for the differences in friction. 


^https://github.com/carlosgs/GNBot 
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